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ABSTRACT 


This thesis considers two models for the computation of 
position finding confidence, one of which utilizes a bivariate 
normal distribution of region, and the other a chi-square 
distribution. The two models are based on different 
assumptions, these are explained and explored. A computer 
Simulation model is presented which utilizes both position 


finding models under varying conditions. 
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Le)6UCA NNT RODUCT TCR 


Bearings are taken on an emitter (target) of unknown 
geographic location by n stations of known location. The 
best point estimate (BPE) of the target is then calculated, 
along with a confidence, or search, region within which the 
target may be said to be located with a given probability. 
Meese cOnocitutes the classic radio direction finding position 
Mecation problem. 

This classic problem is part of a general class of 
position estimation problems whose basic approach is identical, 
wewen 1S the calculation of a point location from a number of 
iemes Of position. 

Variant problems in this class occur when lines of position 
are established from ranges, or relative times of arrival of 
an emission, or when the object to be located takes the 
measurements on the reference stations. 

While this paper examines the classic problem, much of the 
methodology applies to this more general class of problems. 

The BPE is the primary result of a poSition location 
solution. The confidence region amplifies this information 
in two ways. The size of the COntiIdence Teg On~esem 7c cco me 
Measure Of the reliability of the BPE, a Jarger regen 
implying a more tenuous location. The region also describes 
an area in which the target will be located with specified 
¢onfidence. The geometric details of EhisS region peevilae se1] 


more Likely locations for search. 





Li... COMPUTATION OF UCONBIDENCE Re G Lois 


A typical [@} sequence OmPOperatious fForva posmeren 
estimation by computer is as follows: A first point esti- 
mate 1S made and am bane Surface tangent to the earth is 
established at that point. Remaining computations are made 
in the plane rather than the surface of the earth. Unusable, 
or wild, bearings are rejected and a BPE made utilizing the 
remaining bearings. A confidence region 1S computed using 
the same bearings. 

More than one mathematical method is available for most 
of the steps [3, 4, 8, 9, 11]. The two methods used herein 
for confidence region computation are designated the 
bivariate normal (BVN) and the chi-square (x7). 

The Bivariate normal region [3, 4] is the most commonly 
used. This region is the Sane of concentric ellipses 

ax* - 2bxy + cy2 = -2 log (1 - P) 


where P = probability that the target is found within the 


ellipse 
x, y = coordinates in local (to the target) reference 

Syiscem 
n Dye 

as cos 8 5 
a o2 

J 

n Sin 6. cos 6. 

as — tt 
j=l aa 
n sin? 05 

c= 2 5 
j= Or 
j=1 ; 





6. = bearing observed by station j 


variance of displacement of 4th line of bearing 


n = number of bearings 


The bearings from the stations are assumed to be normally 
distributed with a mean of zero, and are displaced parallel 
to the true bearing at the site of the BPE. The confidence 
region lies ina plane tangent to the surface of the earth 
at the BPE as calculated by the least squares method [3, 4]. 

The BVN confidence region requires several assumptions: 

(1) The position lines are great circles on the 
SULrtdaCesor seine wearcEl. 

(2) Measurement error is Gaussian with mean zero. 

(3) Dasplacement of position lines from theseuie 
position line iS parallel to the true position 
line. 

(4) The earth is flat in the vicinity of the 
confidence region. 

These last two assumptions represent approximations that 
are made to Simplify the calculation of the confidence region. 
Thus, the calculated region iS an approximation. The parallel 
bearing displacement is an approximation of the actual dis- 
placement of a point on the bearing by a bearing error. From 
this, the probability distribution of the points is developed 
[4]. Placing the resultant curves in a plane adds an 
additional i Stemson which grows with the size of the 


confidence region [7]. 





The chi-square confidence region [3, LO} 2s"the fan 


piECOncentric curves 


x2 = . i (0. - 8 yo 
j=l $4 : 
where 
s = bearing variance of station j 
. Sst = bearing observed by station j 
B = true bearing from station j to a point 
Cay) in the curve 
n = number of bearings 
x* has a chi Square distribution with n degrees of 


freedom. Of the above listed assumptions reguired by the 


BVN method, only (1) and (2) are required for the Mia 


region. 
The region lies on the surface of the earth centered at the 
minimum x? point, rather than in a plane. 

Both methods require the estimated variance of the 
bearings to be used in the calculation. The effect of an 
incorrect estimation could be much larger than the errors 
See approximation. Estimates in practice are a function of a 
Seaeions past performance .plus the direction finder operator 5 
evaluation of the reliability of each bearing he observes and 
me Ports. 

The assumption of Gaussian bearing distribution eliminates 


the case of wild, or outlying bearings which occurs  praceme., 


largely through operator error. As a Gaussian distribution 





is basic to both the BVN and x2 methods, elimination of wild 
bearings from the computation is necessary, though the theory 
of each method would be affected by their assumed existence. 
Wild bearing rejection methods are available [6, 7]. The use 
of such methods is usually limited to four or more bearing 
Situations. With two bearings, no basis for detection of 
Wild bearings exists. For most three bearing situations, when 
the existence of at least one wild bearing is indicated, it 
is impossible to identify it, or them. With four bearings, 
one wild bearing can be seen to miss the BPE area formed by 
the other three. 

iio pracerce (6, 7iMeehe Outlier detectionwand Trejeceron 
method is used, but a Gaussian distribution is still assumed 
for those remaining. 

Daniels and Vajda [3] object to the use of 2 regions on 
Meeegrounds that it may occur that no points can be found to 
Satisfy the resultant equations. 

Read [10] examines the simplifications required in 
calculating the BVN regions, and compares the two methods. 
The x2 regions are expected to be too small, and the BVN too 
large, though the full effect of the simplifications by Kukes 
and Starik [4] is not known. 

Comparison of the two methods by computer plotting 
techniques shows considerable difference in size. Additionally, 
the effect of the approximations used in the BVN method is 
apparent. While BVN region is elliptical in all cases, the 
x? region may assume an egg shape. 


8 





A barrier to the use of 2 


regions is the lack of a way 
of reporting them. A BVN region can be described by its 
center, angular orientation, and length of semi axes. No 
such simplifying parameters exist for y2 regions. One 
possible approach would be to attempt a conversion into a 
partially equivalent rectangle, Square, or circle as is often 
done for BVN regions [l, 6, 7]. Since these regions are 
always larger than the computed region, contain areas of 
lesser probability, and omit areas of greater probability, 
any advantage of using one method over the other would 
probably be lost thereby. 

A thorough investigation of the two methods is required, 
in order to determine the actual probabilities represented 
by each type of confidence region. The differences indicate 


that at least one of the methods has error. A computer 


simulation was written to compare them. 





Tit. COMPUTER SIMULATION MODEL 


The appended computer program provides a means for 
Simulation of the two confidence region methods for identical 
maout data. 

Inputs to the program are station latitude, longitude, 
and standard deviation of bearing distribution; target 
latitude ‘and longitude. 

The true bearings and distance for each station and 
target are computed and a bearing error generated and added. 
ims portion contains parts of the Bergh program in Pope [8] 
Poets a program used by LT Paul Winkler, modified here to 
include more general positions. 

A first-point estimate of the position 1S made by Pope's 
(9]} vector method, and the final fix by the least squares 
method of Kukes and Starik [4] in a plane tangent to the 
meer at the first-point estimate. This portion was used by 
Lunde [5] but contained errors since corrected. 

A bivariate normal confidence region is computed by the 
method of Kukes and Starik, and orientation angle and 
length of major and minor simi axes generated. Additionally, 
@meput by graphic plot is available. 

A x2 confidence region is then computed and a graphic 
Peet generated. No paramneteric OUtpuc Lor sene x? region is 
available, as no method yet exists for describing the shape 
®f the region. This portion was originally weiecen oye 


Winkler, and since modified to include more general positions. 
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The subroutines GAMA, for the normalized incomplete gamma 
function, and MTRMAP, the graphic plot, are Erom the library 
of the W. R. Church Computer Center at the Naval Postgraduate 
Benool. 

The graphic plots are on a grid measured in degrees of a 
great circle, resembling a Mercator projection at the equator. 
Hence, visual distortion of the regions occurs as their 
location approaches the poles, but the geographical points 
of the ia regions are correct. The geographical points of 
the BVN regions are correct, except for the distortion caused 


by the flat earth assumption. 


ial 





lV... SUGGESTED sSIMUDATI ONE LOer bun 


An undetected error in previous work invalidated all 
data generated. Time constraints prevented regeneration of 
the data after correction of the error. 

The behavior of the regions generated by the two methods 
differs with the number of stations in the net, configuration 
of the net, standard deviation of the error used, and relative 
itecation SE target with respect to the net. Thus, a general 
case may not exist. A valid comparison would include 
sufficient combinations to detect all possible patterns of 
behavior. 

A useful measure of effectiveness for the regions would 
be the difference in the probability level in each of the two 
regions occupied by the true position of the target. 

The standard deviation used in generating bearing error 
is also used in computing BPE and confidence regions. A 
minor modification to the program would permit the examination 
@e the consequences of incorrectly estimating the actual 
bearing distribution of a station, by computing these with 
a different standard deviation that that used to generate 
me Crror. 

The errors generated are from a normal distribution, and 
include no wild bearing. A means of detection for wild 
bearings is not included in the program. Such a routine 
may be added to the program, along with a qenerarien reugln— 


mom wild bearings. 


diez 





The modifications above would permit examination of more 
realistic situations through the use of operational data, 


along with a model of the actual net which supplied the data. 
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V. AN ALTERNATE APPROACH 


An underlying assumption in the methods discussed icetaae 
memere 1S a unlLiorm prior probabilacy for Eh= targee locacr1on. 
That is, the mathematical solution of the BPE and confidence 
regions are inputs to the solution of the operational 
question of target location. Prior knowledge of ‘probable 
meeeeget location does not affect the solution. This approach 
1s intuitively correct generally. 

Operationally, there are two uses of position finding 
miiformation : 

(1) To determine the location of a target of known 
Sr unknown identity; 

(2) To identify a target from a set of known targets 
bveetes location. 

This second use introduces prior information into the 
problem, and occurs usually when attempting to identify a 
fixed land-based station. 

Buecerly {2] describes a method for incorporation of 
information of this type in producing a Bayesian position 
estimate, though he implies that its use could be more general. 

A further modification of the program would permit a test 
of the efficacy of this method by comparing the number of 
correct identifications when using the Bayesian method with 
the number when selecting the target in the lowest level of 


2 


the confidence region. Both the BVN and y metnoes coulame- 


wsed in Gach. 
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MAAONNM 


SOOO? RIT 


AAM 


READ(5,100) NSTA 
100 FORMAT(I2) 
DO 2 I=1;NSTA : 
2 READ (5:102) STALAT(1), STALCN(1}, STDEV(1) 
102 FORMAT(3F10.4) 


READ TARGET POINTS 
ae eet OOF NTAR 
DO =e TAR 
1 RPEAD(S « boo) TARLAT(I) + TARLON(I) 
READ GRAPH SIZE 


READ(5+102) CGRAPH 


PART TWO — GENERATICN OF BEARINGS AND BEARING VAR DANGEs 
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ive COMPUTES BEARING ANGLES AND BISTANCES 
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CHECK FOR SAME MERIDIAN CASE 
meenA.Gl.Js91) GO TO 105 


BB=0.0 

MmeCSTALAT(1).GT.~TARLAT{(J)) BB=180.0 
THETA(1I.J) =BB 

DIST (1leJ )=ABS({ STALAT(I }-TARLAT( J) ) *RADIUS#RADDEG 
eeei0 23 


moe CONTINUE 
CC=ABS{( AA-180. 0) 


16 





ae 
ees 
—_ 


oe 
— 
Y 
© 
5) 
kis 
© 
-— = 
us < 
Oppo.) fod 
Liu Uy LW Ut <{[ C5 
OGe -— = 
i 5 — 
qIet<{< Fe ee) Y 
~ Were. XO oad = 
oud KH HK 4 asec 
lw > LV LVN WY) z = 
3) (ey) eects e © uw 
= WwW Ww Oc Fe Ie a = O 
Oo De ee) —+ <dt> as tty =a} es ane es 
Oo —_O 2) O Wie O@e' = eo Nk yen) <{ <{ 
— FO <{ = ten SSS! & +4 Fey co -— — 
Ie fod qq Nnszweww ade Ce OOo wom Y 
O sil —+ — Chip we YY at ee Of © 
re CO FM Oo tN OamMaad a Om ++i! = — 
SO) a FEF mnjt uw tu 9.09 =) ‘®) 
Or J vee aed 0 Se am oe mom FO fad me 
co pom 2.0) Lic elU<ct <f NNO =) ii tl i) CQ Je eo 
~O Det OQ ae DRR Wu NZ <_ Ge a ae mn. 
—~ mee NY wt J OY) acl Jad ol Oe oe = Nk ko a2) -— 
4 =O —+ ge ett | am co aa wes OQ On 2 <{ 
OD F ONO SS. => | <{<{<I <{ a =) aan ot eM OR, k- 
en ee | (ea a oe 1 nk 8 Pa oP all oe = ree @& eO jl << 
Oo <{ em] i =) <CSCjeirs NZ Z QH ~wt © ~~ OS Onia tu as 
PE oe ee? 1 Lams ma ODDO ODORAR OW NANI SF ng ee e.-@- 9a Cc) O 
re Nege eamMUu O no mw OOnNnM tH xX ae me OOKkkF Ukr e ad 
OO WgoanoenND WM iT ean <ictk e a). Sea u. 
COhea~ONn 2 ty JI 1 Hou Wea kt st F< iy <I © ef & Ome TW 
OeqData~On W D O O00 aAau oJ ea) SS) Se) We) 
QOINE WERE D JRA>PW ANNAN VO Iti OO NRNRRRRke fF OO 
wa SR WON 2 2 ake NMS] Wow —Inmn me VA mee J et er LL iI ig 
LODLTWRHOO ws >POwrwm OORN Am FFO << ed Se as eae <{ 
MOV SNOW Fe BOG DOA o>) a Wo NR RRR OF - 
Fd ow lw WY 
Oo © © 6a a 
Oe op) -— CO 
| 
OUW O OO S) () fad 6.9.6) OQUOUUOO0 O 


GAMMA(I,J)*RADDEG 
my 


GAMMA( 1 .J) 





STALAT(1) )+COS( TARLAT( J) Y*CUS(SIARAI 


DEST ANCE iN NAUTICAL MILES 


23 CONTINUE 


LST 


C 
a 
3 


WO UOUUUUOUOO 


ee) 
Oo 
Ww 
~ — 
© @(5 
2 HL 
~ (ae) 
+O 
+ pag 
tefee 
Ce x % 
Kw OO 
NSN Al 
= ‘> a ON | a 
~ (OG) 
ae e e<t 
(a a 
Nt wo | t x 
ON NO ia 
e «+O 208) wy CW) 
Oe a 0 © Je grea 
Saito “NiO (C Oe 
aeQ eM & EM 
NOrt Se Ss oe 
ae NO oe il 
—> ™ CJ NY OO 
rn ae os. > ap | “— 
> Ow Om @ e (ne) 
tlim—w © eQOal vm ee 


QOhMSR AS ee ALN Qiawr US 
KFmNLOMATRH OOM YY ~D 
NeITn [ia A Ae-= 
Oe Occ CHRO 
boom hem HE HR He 
I Be ree 
A. trOOiwLLOooooO GO; 
NZL NORRNONIDNNO 


N OFO DB Odo 
N ANN OM NO 


a, Ee 


1+NSTA) 


36X%+45H BEARING ANGLE FROM ITH STATION TO Sia. A RG Eales 
L1sNTAR) I 


AM MST A One tO sit A CG tallesy a7 
NTAR) + IT=1sNSTA) 


ls 


cel 


18 





ey Pe 
EL se eS 
f-t- cal 
<[ Ge 
a=) <{ 
gat oe 
Zt (=) ale 
m4 <f \0 = 
Y) ~~ -— 
eS <x 7 
<= ao~ 
<[ ae ~~ 
SS © a Cc m4 
faa bh Ts Wa 
réO = 
oe) — << ae 
<I 0 O 
‘a dta's <I — 
Sn) b— onl - 
(8 igre xc“ MFR uJ 
Lich e jm <I 
> a Ga ze 
Gaal 6 Oe ae amm 
eo (0 = oY il I~ 
— 7) ~ e <TR Y) 
CaS, _ = ik uJ 
<1 => o Gy - © 
eal Y) ae Te ee 
a Uy ©) ag 4 i on tO es, 
a - UJ © — a — 
GO < Om ‘ad Ya 7 — 
EU a) — OrX< Oe age GO pre — 
= [ee <I eo WJ Otsu © Y 
vat << (OB aS, (3<I rH ©) 
(ah a > “LU an) eh Il Ge 
= -— © LJ ~—~WWY7D 
O<ct mad “~ - (© OP BRS Le 
Se e <I a) <I LU\——~ & © 
oO =. — © >< —| fon mon 
Ges iad wey “fp ) o> Pad 
a C3 bon © <{><hk- e-= © 
eae 4 a Ne 7 Y) aoe Eee, 4 Ge _ 
<1 ale LLY) Y) <I Jun ke 
k- —sstL uy IU OK << t «<{-- = << (as 
wwzma 4 Fe LL o Make MU) -— <_ 
IO a Yr ale ae, EE =) b— 
oO ed < Iver Ne o oe OO em UN a Pea 
t~ Oo Reet ei oe — oH Hd = ° 
ILD =) Sole eee Ss Tem YU. © 7 
OMe De Geert k YO Oo<tO o~ Q Hi 
= on o oe )~ jl a) YM ONFWUAO al 
ZOW NN aa ene NN tc a AI Wise | i) 
Or ied Ween 7 ee LY ae) owe Te oN 
mUAb- OO mM ew CO D 5 O KRORM~ us rey 
aed al —br rOn =a 2 a eo LJ N 
WOOO Wat @o@oo) — a qj wrtowq< Ge an) 
Uy ke = srw t ke F&F qe en! te ag 2 
ple | all 6 Hitt =) 2 b— mee tO ae 
AY COO SR aaa) eC aca SD eS) 
NEO Bu ae oe Q Swot u _ ee 
— O Oo |= mod (es 
ae Cr © © eS, oS Ov <{ 
apes ol Oo ey Ge) 2 uy a 
N N ON 


OOOO O ©) Om O OUOUOUOUOU O 


<I 
Ud 
co 


LOCA 


_ 
— 
_— 


KEN ON A TARGET 


I= 


BETACIeJ1),THETA(I,J1)+ 


i 
12XsGHLONGI TUDE »L4Xe5SHTHETALIXyL2HTRU 


I) ,STALON({ I), 


aOUY om 


COWYOZzrtO 
TS LL be SLL 
| 


iS 


ep) 
MN 
N 


Te, 





FIRST POUNT ESA. 


A 
E 
POmPUTE BEARING PLANE INTERSECTIONS 


INT ERS CORTONS., 


LINE 
ITTUDE s4Xy 


WMOMNOIMNMO 
| 
Te | 


PS Ts 
OO 


ae ee 
On 
3% 
KE 
Ig< 
_ ae 
OUWN 
Mow at 


LL LY tJ 
>< NM 


SBNG* SLAT*CLON 


mC oor CMCGN = SENG@SEAT=S EON 


CBNG=SLON - 
SENG 6 LAT 


a im, 
metho Pt 
eet See? eet 
Loud 
<OO 


LJ UL tL) 
> NOX 
ie K 
_= —— 
— rt pn 
ett Sot eel 
CE TE ie] 
Ota 


Lust) 
NIK > 
eM 
= 
Pt rt Pe 
ee nati ad 
Lig Li Lu 
MOo< 


Hoo a 


— a, 
— | Pe 
wae et Set 
LW uy UL) 
Ou tL 


| 


88 CONTINUE 


OOO 


OOO 
i | 


=a Ee 
===> 
NYY 
= 


J 


es 
LefLIMIT 


Jz 


Last! 


OO aio 
H 36 4 
ms oO wae 
rt peg 
ent et eet 
tL) Wu) 
MDO<cxt 


LU LU LL 
D> = 


20 





GO 1G 9o 
GO TO 91 


AO 


WLW LY 
a> = 
ito 


LU WL Li 
De 


@eMPUTE CENTROIO OF INTERSECTIONS 


WaLSSTG) 
a> = 


++ + 


paella 
>=) 
MWY 
I> = 


Woe 


ee al 
ee Fan 
MMV) 
D>= 


ee 
ON 


1293%eF 84.394X9F 8. 3} 


LCN 


TUE A) a 
Loa) 5 
BEARING PLANE INTERSECTIONS 


J 
5: 


ULON 


ULAT 


WOO 


Te 


eS .es 
M+VSUM*EVSUM+4W SUM*WSUM))/RADDEG 


g 
id 
Zl 


ASUM 


CONVERT CENTROID TO LATITUDE AND LCNGITUDE OF FIRST POINT 
Baum 


BO CONTINUE 
A INUE 


89 CON 





a) 

ii} 

QC 

cS O 

ang LLy 

(we ee 

+ QO 

— << 
It <ae 
rw FH 
ANS i 


CO<tt- Ov ii I 
woz tte 
OoddMNMYM 
ONUWO «KX 


8999 


IN LOCAL REFERENCE PLANE USING NAPIER'S ANALOGIES 


VY) 

z= 

<I 

rt 

Ee) 

<I 

Cc 

Tae 

— 

O 

Fae 

imi 

oe WY 

<< (ous 

uJ pane 

wm bet 
e ot) 

Ze <{t 

C) LEC 

m coe 

ke 

<[ ij = 

I E<t 

Y) a 

il Oo 

i C2) 

VY) @z. 

~< Ud 


OOOU0O VO 


)*RADDEG 


ROS 2. OFAN a Rin’ 
NORTH OR SOUTH BEARING 


CHECK FOR DUE 


OVUWO 


ce 
® 
© 
th 
4 
— @ 
LN ~ 
oe ee (UU Oo '@) 
mM NS or il 
an) ont 
Cc) om — 
a o ~ 
<t tH << 
oO -—! <I 
CP er, =e 
~< OO AN 
~ e be 
© e mH re *® 
e ~~ (A>) we C) 
© — jee r- @ 
fs Gufe= be Oe 
Ov efittack li<f esr 
em MDM Oe LUM 
Org 2WNY tin 
Tee OMMY~ lOO 
COIN oativuakr 
Sew ZI I~ 
Hot, Cea thw) oe 
mI OO INO 
uy 
vy 
~o 


om) 7) ~ 
= ee N 
wa © - 
ge, 1 GC) = =— 
<< > <a 
ei CL Ww vee 
oe —~ w 
E- = 2 
a << oO 
ib reOM ho ~ 
2a ae % De 
=e mon OL 
Zz Sea A> <t | 
QO Ire HOw 
a fw =2ANo 
Ik >< OK ss DLW a” 
DQ New WOW 
SS ie SO 
Lo ae tA oa 
> WW O8ee ss eet 
=< Ov db NMA 
— tly) SST MAIN Hiatt 
= 2 Steaia = <i I 
re el em Pa 
O a Wie) Awe 
OO KF ygarrt aAoae 
Fras 
Oo O 
fF OO 
fad) 
OO OUO W 


ZZ 





TANCE FROM LOCALIZED BEARING TO 


Fes 
OQ 
oa] 
ke 
Prd 
us 
> 
Pog 
© 
SS 
z 
5 
ad 
VY 
5 
ae 
ad 
t 
a 
-_ <{ 
i = 
_ ce 
<I © 
- = 
Li Tp) 
NI @ Zz 
| a ee, 
iN ra 
ad = O r- 
a te Ge tu <{ 
SS Ov LL) se Lid a) = 
U\— Ne) ~ o © Ov Ue<t > oe) mB 
LAN uN N Oo e Or UF Pad <I 2) 
a8 Kona) () a OQ Oo Ov NW iu ©) Ys Tes ~ al 
= = _ e fi on) aN ) pas LL ae ~ <{ 
ee) G) =) OO ea Oy — Ok - —4 S 
— ~~ to ete <t Oo IO j— Ye > LL <f C5 
—Oz OC Oa CO mn<z 1 - @ ap. <V) LL z= 
NO << ap) <{ © AIX OS LL Ls © Pa <f pee | <{ UJ — 
“SF - =Z2 lo hee Nie - = eS <IN = 
ae i ~ me DRaAZR Lis - a Lis os Oe i Be O 
fee sa = am ZHAO a a ae <{ mln mf ff Ta =e 
Lon ors sz Oe DOM"O —~t. WW e Liu —_ O>7 sera —J 
Or Ome) esq) waIOM}t *=— OU OW cO Oa ze el) @ 
+ o> ne 2 =a Nt enmOow a 6 OO O fF iO eso. tL 
P<) Ont DW ZAOW lwWwe~w | en ej Ua auyee |<t>ke 
wD OS ET ae, es YO az ewe asa St = Lj WwW MN ‘oa | ect) f: 
CRE ewORFO~M Se FRRFOIKWWr. Oe N > | wo Ow mm (1 
HQNONMVONNWNANW IYI me NOM Il om 10 on eit a0 7 © Le 
MeO *M @ COMDOODINE eZWZ0 eo mnt DD eEe- ~—~0O ~~ 
onto Faas ZOUZN SHH HOU <f salt oc <f <I we WL Nase oe) 
Ire hv hOrFNNOmUHReZi MNMNEtNOw tl wells UE), oO ~— reO Zilkeo =e 
WOM OR CE ORFS OrREUX ywonnubu Ato bf Bi) Ly << <I rat LL) OD ems ‘aa 
QOw32 wDew 2522) 28 See ee | Zz Fe JY Cui WwW 
RPL ZWwOLZwoOZwOor¢troumuLnro Lie. jes WwW >W eto = 
YCHMOORnMNOOUNMODS aIaTReOnWihL reo Nees O << Ww ™N Oe wea 
- < ae > 
uy O Mh WwW < <{W =) 
Ww Tal Tal +s-ON tr OO Y 
faa) an) an) ay 


O QO OQOUMODO  QOOOOUL OVUOO 


78 





CALCULATE ANGLE BETWEEN MAJOR AXIS OF ELLIPSE AND MERIDIAN 


ASUM=A SUM+ (COS(ZETACI) )¥COS(ZETA(I)) )/(SDRAD(CTI)*SERAD(T) } 
BSUM=BSUM4(SINCUZETACI) )¥COS(ZETACI)) I /(CSDRADC(I)*SORAD(I)) 
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a CONTINUE 
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ESUM=ESUM4+Pu (K)*( (ASUM* SIN(ZETA(K) )-BSUM*COS(ZETA(K) 1) 7 (SDRAD(K) 
Q L1SDRAD(K) 2) 
°° CONTINUE 
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4 Wiese =X1LAT+YO/ RADDEG 
G Mesde YEST = LEAST SQUARES ESTIMATE OF POSITION 
o 
WRITE(6.9007) YEST.XEST 
MOT FORMAT(* *.'THE LEAST SQUARES ESTIMATE OF POSITION IS ° yF6.2s 
4 DEGREES LATITUDE'//45X%eFf.2¢' DEGREES LONGIIUDES///} 
9 
C 
; 
. PART FOUR — COMPUTATION OF CONFIDENCE REGIONS 
4 
Mme THIS SECTION OF PART FOUR CALCULATES A BIVARIATE NORMAL CONFIDENCE 
5 REGION (DANIELS, KUKES AND STARIK) | 
. 
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N=13 
M = 13 
ey = 0.0 
Be = 9.0 
BND = 0.0 
PMIN = 0.0 
TJT = 0 
ICON = Il 
hex = i 
DO 59394 II = 1424 
‘ Boe T{(J1) = 0.0 
meer THE Q LEVELS 
CALL MT RMAP CY, NM, Ts BND AZ »,BZ,AMINeISJT,ICON,IGR) 
WRITE(6,592) 
Bo2 FORMAT(///-.45X. "LEVELS OF CONSTANT VALUE OF Q*) 
: WRITE( 62509) CGRAPH-eCGRAPH 
PLOT THE CHI-SQUARE REGION 
CALL MTRMAP(CeNeMeT»eOND AZ eOZeAMINIJT,ICON,IGR) 
WRITE(6.593) _ 
Mee PURMAT( ///+45X -'CH1-SQUARE CONFIDENCE REGIGN! ) 
| WRITE(6.509) CGRAPHsCGRAPH 
<! CONTINUE 
STOP 
END 
BeesOULTINeE GASS(CIX.sSeAMeV) 
P= 0.0 
Bie 50 L=1-12 
Grek WRANDUCIX, 1Y.Y) 
Ix = 1lY 
mo A = AtyY 
V = (A~6.0)=St+AM 
RETURN 
END 
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